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Abstract This paper presents a new theory of the dynamical tides of celestial bodies. It is founded on 
a Newtonian creep instead of the classical delaying approach of the standard viscoelastic theories and 
the results of the theory derive mainly from the solution of a non-homogeneous ordinary differential 
equation. Lags appear in the solution but as quantities determined from the solution of the equation 
and are not arbitrary external quantities plugged in an elastic model. The resulting lags of the tide com- 
ponents are increasing functions of their frequencies (as in Darwin's theory), but not small quantities. 
The amplitudes of the tide components depend on the viscosity of the body and on their frequencies; 
they are not constants. The resulting stationary rotations (often called pseudo-synchronous) have an 
excess velocity roughly proportional to 6ne 2 /(x 2 + x~ 2 ) (x is the mean-motion in units of one critical 
frequency - the relaxation factor - inversely proportional to the viscosity) instead of the exact 6ne 2 
of standard theories. The dissipation in the pseudo-synchronous solution is inversely proportional to 
(x + X -1 ); thus, in the inviscid limit, it is roughly proportional to the frequency (as in standard the- 
ories), but that behavior is inverted when the viscosity is high and the tide frequency larger than the 
critical frequency. For free rotating bodies, the dissipation is given by the same law, but now x is the 
frequency of the semi-diurnal tide in units of the critical frequency. This approach fails, however, to 
reproduce the actual tidal lags on Earth. In this case, to reconcile theory and observations, we need to 
assume the existence of an elastic tide superposed to the creeping tide. The theory is applied to several 
Solar System and extrasolar bodies and currently available data are used to estimate the relaxation 
factor 7 (i.e. the critical frequency) of these bodies. 



1 Introduction 

During the 20 th century, many versions of the Darwin theory, or of what has been called Darwin 
theory, were used in the study of the tidal evolution of satellites and planets (see reviews in Ogilvie 
and Lin, 2004, Efroimsky and Williams, 2009). Those versions were not exempt of problems, the more 
important appearing when they are used to study solutions near the spin-orbit resonance. All of them, 
consistently, show the existence of a stationary solution, which is synchronous when the two bodies 
move in circular orbits or super-synchronous otherwise (that is, a solution in which the time average of 
the rotation angular velocity is constant and slightly higher than the orbital mean motion; it is often 
called pseudo-synchronous). The excess of angular velocity in the stationary solution when the bodies 
are in elliptical orbits is physically expected. The torques acting on the bodies are inversely proportional 
to a great power of the distance and, therefore, much larger when the body is at the pericenter of the 
relative orbit than in other parts of the orbit. As a consequence, the angular velocity near the pericenter 
will enter in the time averages with a larger weight and will dominate the result leading to averages 
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larger than the mean motion n. Then, in the case of a planet or satellite moving around its primary in 
an elliptic orbit, we should not expect the synchronization of the two motions. Pure tidal theories lead 
to synchronous stationary solutions only in the circular approximation. In standard Darwin theorie^l, 
in the usual approximation in which only the main zonal harmonic is considered, the average excess of 
angular velocity of the body is given by ~ 6ne 2 (e is the orbit eccentricity) (See Goldreich and Peale, 
1966, eqn. 24). This is a quantity independent of the nature of the body and is thus one important 
difficulty of these theories. 

This prediction is, however, not confirmed by the observation of planetary satellites. Titan, for 
example, should then have a synodic rotation period of about 8.5 years (i.e. ~ 43° per year) while 
the radar observations done with the space probe Cassini over several years showed that the actual 
rotation differs from the synchronous spin by a shift of ~ 0.12° per year in apparent longitude (Stiles 
et al. 2008, 20100. In the case of Europa, the value predicted by the standard theory is less than 
20 years. The present position of some cycloidal cracks confirms a non-synchronous rotation, but the 
comparison of Voyager and Galileo images indicate a synodic period less than 12,000 years (Hoppa et 
al. 1999, Greenberg et al. 2002, Hurford et al. 2007). The most striking case is the Moon whose average 
rotation and orbit are synchronous notwithstanding an orbital eccentricity 0.055. 

The only way to conciliate theory and observation is to assume that in all these cases an extra 
torque able to counterbalance the tidal torque is acting on the body, and the most obvious assump- 
tion concerning this extra torque is the existence of a permanent equatorial asymmetry of the body 
(Greenberg and Weidenschilling, 1984). This can explain the case of the Moon and a quick calcula- 
tion including the tides raised by the Earth and the C31 component of the lunar potential result in 
a spin-orbit synchronous solution; the net effect of the tides is just a deviation of the symmetry axis 
which is not pointing to the Earth but shows a small offset (see Ferraz-Mello et al (2008), eqn 46). 
In the case of Titan and Europa, the situation is more complex. In both cases, we may assume a 
permanent equatorial asymmetry, but the resulting solution is then an exact synchronization, not a 
slightly non-synchronous rotation as some observations seem to show. 

Several attempts were made to explain these differences either by assuming ad hoc mass distri- 
butions inside these bodies, or by modifying the Darwinian theories. Since the 6ne 2 law results from 
the theory when the lags are assumed to be small and proportional to the frequency of the tide com- 
ponents, one immediate idea is to substitute the linear dependence by a power law (see Sears et al. 
1993). However, without physical grounds to support the assumption, the result will depend on the 
ad hoc fixed powers and remain only speculative. Efroimsky and Lainey (2007) have proposed to sub- 
stitute the linear law by an inverse power law. The grounds for their proposal are some laboratory 
measurements and also the determination of the dissipation affecting the Earth's seismic waves at 
different frequencies. An inverse power law brings with it an additional difficulty because any quantity 
inversely proportional to the frequency tends to infinity when the frequency goes to zero. Efroimsky 
and Williams (2009) and Efroimsky (2012) claim that this difficulty can be circumvented, but it is 
done at the price of an extremely complex modeling. The results of this paper (Sec. 1101 1) may help 
the understanding of what happens in the immediate neighborhood of the frequency zero and why no 
actual singularity exists. 

A different approach was presented in Ferraz-Mello et al (2008) (hereafter FRH) in which the lags 
remain proportional to the frequencies, but the non-instantaneous response of the body to the tidal 
potential is taken into account. In that approach, the resulting excess of angular velocity is given by 
~ Qne 2 {k\ / ko) where kg and k\ denote the response factors of the body to the semi-diurnal and the 
monthly components of the tide, respectively These response factors are not equal. In the stationary 
condition, the frequency of the semi-diurnal component approaches and fco approaches its maximum, 
the fluid Love number kf. On its turn, the response of the monthly component of the tide will depend 
on the viscosity of the body. If the viscosity is small, the body will respond faster and fei/feo ~ 1. 
The result is again ~ 6ne 2 . However, if the viscosity is large, the deformation of the body does not 

1 I will have to refer often in this paper to theories derived from Darwin's theory. To use a simple label, I will 
use the word "standard Darwin theories", or, for short, "standard theories", to denote all Darwin-like theories 
in which an elastic tide is delayed by lags assumed small and proportional to frequencies. 

2 A re-analysis of the data by Meriggiola and less (Meriggiola, 2012) has not showed discrepancy from a 
synchronous motion larger than 0.02° per year 

3 We use in this paper the same tide component names used for fast rotating bodies (Type I of FRH) 
regardless of its actual rotation speed. In next sections we use the names monthly and/or annual for tidal 
components with the same period as the orbital motion. 
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attain its maximum theoretical extent; that is, ki/ko < 1 and the resulting excess of angular velocity 
is smaller than 6ne 2 . However, as in the discussion above, without physical grounds to support the 
chosen value of fci, the result will depend on ad hoc fixed values. 

A difference in the response factors was also considered by Darwin (1880), but it was only spo- 
radically been considered in some papers (e.g. Alexander, 1973; Wahr, 1981; Efroimsky and Williams, 
2009). However, the differences among their response factors were not sufficient to solve the prob- 
lem highlighted above, mainly because those differences disappear when the so-called "weak friction 
approximation" (Alexander, 1973) is introduced. 

One structural difficulty comes from the use of Love's theory of elasticity. The use of Love's theorem 
as a shortcut to obtain the potential of the field spanned by the tidally deformed body without having to 
calculate beforehand its figure of equilibrium introduces the constant Love numbers as response factors 
for all terms issued from the same spherical harmonic of the tidal potential. The only free parameters 
are, then, the ad hoc introduced phase lags. After Darwin (1879), the tangents of the phase lags are 
proportional to the frequencies of the tide components, one result also found in this paper. However, 
since the lags are introduced in standard theories as arbitrary ad hoc quantities, different laws can 
be postulated. They can be freely postulated or, as in some more elaborated investigations, suggested 
from the study of delays in damped oscillators (see, for instance, Greenberg, 2010). 

However, in most of the standard theories, the viscosity of the body is never explicitly considered 
and strictly speaking the so-called "viscoelastic" approaches are more adequately described as "elastic" 
and "delayed" , their actual ingredients. 

We did not consider in this short account theories based on energy dissipation instead of phase 
delays, because they follow a different line of thought (for a comprehensive account of them, see 
Eggleton, 2006, chap. 4, and Migaszewkski, 2012. See also Bambusi and Haus, 2012). The results issued 
from these theories are formally equivalent to Darwin's theory when lags are assumed proportional 
to the frequencies and lead to the same pseudo-synchronous results discussed in the beginning of this 
Introduction. 

This paper introduces a new rheophysical approach in which the body tends always to creep towards 
the equilibrium by the only action of the gravitational forces acting on it (self-gravitation and tidal 
potential) and does it with a rate inversely proportional to its viscosity. The adopted creep law is 
Newtonian (linear), and at every instant the stress is assumed to be proportional to the distance from 
the equilibrium. This leads to eq. ^ and with only one exception every result in the paper is a direct 
consequence of this first-order differential equation. The exception occurs when studying the shape of 
the tide in hard bodies as the Earth. In this case, in order to conciliate theory and observation, it is 
necessary to assume that a purely elastic tide (See Sec. [TU]) exists, superposed to the creep tide. 

The physical model is presented in Section 2 and developed in sections 3 and 4. Section 2 also 
presents a short application to the case of bodies in circular motion with results equal to those obtained 
by Darwin (1879) and which served as the basis for the introduction of a lag proportional to the tide 
frequency in his 1880 theory. The next sections are devoted to the calculation of the perturbations: 
First, the perturbations on the rotation of the tidally deformed body (Sec. 5) and its synchronization (or 
pseudo-synchronization) (Sec. 6); then, the perturbations on the semi- major axis (energy dissipation; 
Sec. 7) and on the eccentricity (Sec. 8). Section 9 discusses the value of the relaxation factor 7 on the 
basis of our current knowledge of the tidal evolution of stars, planets and planetary satellites. The 
creeping tide theory is completed, in Section 10, by the introduction of an additional elastic tide, 
necessary to reproduce the observed shape of tidally deformed bodies. 



2 A simple rheophysical model 

The usual standard model considers an elastic tide and delays the tidal bulge by an ad hoc phase lag 
in order to take into account the body anelasticity. In this theory, we propose instead of that, a simple 
rheophysical model, the basis of which is shown in fig. 1. We consider one body of mass m and assume 
that, at a given time t, the surface of the body is a function £ = (({p*,9*,t) where £ is the distance 
of the surface points to the center of gravity of the body and (p* , 0* their longitudes and co-latitudes 
with respect to a reference system rotating with the body. In the same instant t, the body is under 
the action of a tidal potential due to one second body of mass M situated in its neighborhood. No 
hypothesis is being done on the relative importance of the two bodies. Both may play the role of the 
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Fig. 1 Elements of the model: £ is a section of the surface of the body at the time t; p is a section of the 
surface of the equilibrium spheroid at the same time. 



central body and its satellite or hot planet. In actual applications, both cases have to be considered 
(see FRH Sec. 18). 

Would the body m be inviscid, it would immediately change its shape to the equilibrium configu- 
ration. In the simplest case, the figure of equilibrium of it under the action of the tidal potential is a 
prolate Jeans spheroid p = p((p*,6*,t) (see Chandrasekhar 1969) whose major axis is directed along 
the line joining the centers of gravity of the two bodies. If a e ,b e are the principal equatorial axes of 
the spheroid, its prolateness is 



fie 




(1) 



(Tisserand, 1891) where R e is the mean equatorial radius of m and r the distance from M to m. Terms 
of second order with respect to e p are neglected in this and in the following calculations. 
The adopted model is founded on the law 

C = 7(P"0- (2) 

The basic idea supporting this law is that because of the forces acting on the body (self-gravitation plus 
tide), its surface will tend to the equilibrium spheroid, but not instantaneously, and its instantaneous 
response (measured by £) will be proportional to the radial separation between it and the equilibrium 
spheroid, p — (. Eq. ([2]) is the equation of a Newtonian creep (see Oswald, 2009, chap. 5) where the 
distance to the equilibrium was considered as proportional to the stress. It does not consider inertia 
or azimuthal motions, which may exist and should be considered in further studies. The superposed 
elastic tide whose existence stems from the comparison of the observed shape of the tidal deformations 
and the theory, is not included in the equations of the model because the force due to the elastic tide 
is radial and its torque is zero; it does not affect the rotation or the averaged work and eccentricity. 

The relaxation factor 7 is a radial deformation rate gradient and has dimension T _1 . It is 7 = in 
the case of a solid body and 7 — >• 00 in the case of an inviscid fluid. Between these two extremes, we 
have the viscous bodies, which, under stress relax towards the equilibrium, but not instantaneously. 

We shall mention here the very similar equation used by Darwin in his first paper on the precession 
of a viscous Earth (Darwin 1877jj to define its rate of adjustment to a new form of equilibrium, soon 
extended to the study of tides (Darwin 1879) by means of a law similar to Eq. $2$. However, in the later 
paper, he was rather interested in the ocean tides upon a yielding nucleus and was not satisfied with 
the results that he classed as fallacious. For unless the viscosity [of the Earth] was much larger than 

4 We may paraphrase one of Darwin's statements by just changing the symbols used in it by those shown 
inside brackets: But because of the Earth's viscosity, [Q always tends to approach [p\. The stresses introduced 
in the Earth by the want of coincidence of [Q with [p] vary as [p — C] • -Also the amount of flow of a viscous 
fluid, in a small interval of time, vanes jointly as that interval and the stress. Hence the linear velocity ( on the 
map), with which [Q approaches [p], varies as [p— £]. Let this velocity be py(p — C)]j where [7] depends on the 
viscosity of the Earth, decreasing as the viscosity increases. 
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Fig. 2 Equatorial section of the equilibrium spheroid corresponding to the tide generated by M on m. Angles: 
a is the distance from the generic surface point to the vertex of the spheroid; v is the true anomaly of M; w is 
the angle between the origin meridian of the body and the pericenter of the orbit of M; ip* = zu + v + a. 



that of pitch, the viscous sphere would comport itself sensibly like a perfect fluid, and the ocean tides 
would be quite insignificant. This is perhaps the reason for which he never went beyond the circular 
approximation later used (Darwin 1880) to introduce the tide lag and the tide height in the Earth 
model in the study of the secular changes of the orbit of the Moon and the changing rotation of the 
Earth. 

It is also possible to obtain Eq. @ by integrating a spherical approximation of the Navier-Stokes 
equation of a radial flow across the two surfaces, for very low Reynolds number (Stokes flow), a case in 
which the inertia terms can be neglected and the stress due to the non-equilibrium may be absorbed 
into the pressure terms (see Happel and Brenner, 1980). The boundary conditions are £ = at £ = p. 
The pressure due to the body gravitation is given by the weight of the mass which lies above (or is 
missing beloftQ) the equilibrium surface, that is, — iu(£ — p); the modulus of the pressure gradient is 
the specific weight w. 

This comparison allows us to see that the relaxation factor 7 is related to the viscosity coefficient 
n through 

= wR 3gm 

7 2r; 87ri?V 1 ' 

where g is the gravity at the surface of the body and R is its mean radiufl This equation is important 
because it allows us to estimate the range of possible values of 7 for the celestial bodies to be considered 
in the applications. 



2.1 The creep equation 

The function p may be written as 

p = Re (l + ie p cos2<^ (4) 

(to the order 0(e p )), where \P is the angular distance of one generic point on the surface of the 
equilibrium ellipsoid to the axis of the tidal bulge (that is, the direction of M). 

If we restrict the present study to the case of a "planar" problem, in which the orbital plane of 
the tide generating body M cuts the body m symmetrically (i.e. m is symmetrical with respect to an 

5 This does not mean that a negative mass is being assigned to void spaces; it means just that the forces 
included in the calculation of the equilibrium figure need to be subtracted when the masses creating them are 
no longer there. 

6 Darwin (1879) used a very complete construction of the Navier-Stokes equations and, in his results, the 
numerical factor is 3/38 instead of 3/8. His numerical factor is determined by the spheroidal form of the tidal 
potential, but the intensity of the potential does not appear in the result. So, his result would hold even for an 
infinitesimal tide! 
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equator, and M lies on the same plane as the equator of m), the differential equation of the adopted 
model for the creep becomes 

C + 7C = IP = jR' + ipRetp sin 2 6* cos(2£* - 2m - 2v) (5) 
where 9* is the co-latitude (introduced through cos!? - = cos a sin 9*), R' = R e (l — ^e p cos 2 9*) and 

(p* — m + v + a (6) 

(see fig. 2). 

In the solution of this equation, we have to consider that = Q, angular velocity of rotation of 
the body m, which is assumed to rotate in the same direction as the orbital motion of M. 



2.2 The circular approximation 

For the sake of making clear which are the main consequences of the proposed rheophysical model, 
before considering the full model, we consider first the simple case in which the relative motion of the 
two bodies is circular. In such case, r — a (semi-major axis) and v = £ = nt (mean anomaly). 

The resulting equation is a trivial non-homogeneous first-order differential equation with constant 
coefficients whose solution is 

C = Ce- 7 * + J R' + ^cos(2a-cr ) (7) 

where C is an integration constant. A, <to are undetermined coefficients which may be obtained by 
simple substitution in the differential equation and identification as 

(T = arctan — , (8) 

7 

where e' = e p sin 2 9* and v is the semi-diurnal frequency 



u = 2a = 2Q-2n. (10) 

The integration constant C depends on (p* (the integration was done with respect to t) and may be 
related to the initial surface £ = ((<f* , 9* , 0) through 

C = Co - R' -Acos(2a{0) -a ). (11) 
The solution depends on 9* via its influence on the constants R' and A. 

3 Tidal deformation of the body 

In order to develop the theory, we have to consider the two-body equations and introduce 

r=f~^ (12) 
1 + e cos v 

and 

v = £ + (2e - ^-) sin^ + ^- sin 2^ + sin M + C(e 4 ) (13) 
4 4 12 
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into eq. ([5]). The resulting creep equation is 



,2" ' — ^ 3 



C + 7C = 7i?' + ^^ — — cos(2£* - 2w - 21) + (14) 



e 
2 

£! 

1G 



,2 



(^7 cos(2£* - 2ro - 3^) - cos(2£* - 2w - £^J + y ( - 5 cos(2£* - 2ro - 2£) + 17 cos(2£* - 2ro - 4£) 
( - 123 cos(2^* -2zu- U) + cos(2^* - 2w - i) + ^ cos(2^* - 2vj - U) + ^ cos(2£* - 2w + + 



or 
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JV 



C + 7C^7^+ 157J? ; Sin2r (-)( — ) E ^2, fc (e) 008(2^-2^7+ (A -2K) (15) 



k=-N 



where ./V is the adopted order of approximation of the Fourier series and the i?2,fc(e) are the eccentricity 
functions appearing as coefficients in eq. (TT4"]) . They are some of the Cayley expansions (Cayley, 1861). 
An elementary calculation using simple concepts of Fourier analysis shows that 



E 2 ,k{e) = J- f (^) 3 cos(2v + (k-2)£) dl (16) 



r 2n 

2Wo W 
(see Appendix IA"| . 

The integration of eq. (fT5)) is trivial. If we write ( + j( = F(t), we know that the general solution 

is 

C = e" 11 J F(t)e*dA (17) 

or 

C = Ce-* + R e + R»(0* , t) + 15 ^ R Sin2 °- V £ (18) 



k=-N 



vV + + fen) 2 



where a = — tu — ^ (N.B. a — a = u — £), 



i?" = --i? e cos 2 r e - 7 * / e e 7 *cft. 



* Jt 
and 

/ fcn + v \ 

(Tfc = arctan I I . (19) 

It is worth emphasizing that the Uk are not ad hoc lags plugged by hand, but constants introduced 
during the (exact) integration of the creep equation just to allow us to write the solution in simpler 
form. However the a% play a role similar to the Ek introduced as ad hoc delays in FRFfl. From their 
definition, it is also clear that the o~k are not small quantities as the lags are assumed to be in standard 
Darwin theories. 

The solution of the differential equation can also be written as 

1 ~ - 

C = Ce-^ + R" + -R e sin 2 6* ^ e fc (e) cos(2a + U - a k ), (20) 

k=-N 

where we have introduced 

1 3 . . / M\ ( R e V 

£fe = ~7~E/2,k (c) cos o~k I — — . (21) 

4 v m /W 

£ is formed by the superposition to one sphere of the bulges of several spheroids whose prolatenesses 
are the 



7 The subscripts used for the cr*, are not the same subscripts used in FRH for their homologous lags e fc . 
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Fig. 3 Main component of £(tp*,t). Angles: a is the distance from the generic point to the direction defined 
by M' (fictitious body in uniform motion; M' = M when M is at the pericenter); <jq/2 is the lag of the actual 
ellipsoid £o w.r.t. the sub-M' point (i.e. the axis of the mean equilibrium ellipsoid). The figure corresponds to 
the case v > 0. 



4 The attraction of the tidally deformed body 



In order to proceed, we assume that the potential of m is the sum of the potential of one sphere plus 
the potentials due to various ellipsoid bulges the intersections of which with the equatorial plane are 
given by the sum of terms in eq. (|20p . 

For instance, when sin#* = 1, the main component of C(<5) corresponds to 

R e + -R e e (e)cos(2a - (To), (22) 

which is the equatorial boundary of a spheroid of mean equatorial radius R e and prolateness eo displaced 
of an angle Co/2 with respect to the direction of M' (fictitious body in uniform motion; M' = M when 
M is at pericenter). (See FigJ3]). 
The equation of this spheroid is 

Co = Re f 1 + ~e cos 2#o) = Re (l - ~e cos 2 0* + ie sin 2 0* cos(2a - <r Q ) \ . (23) 

where <?n is the angular distance from the generic point to the vertex of the spheroid. 

(o differs from the component selected in eq. (|2TJ|) by the additional term — cos 2 0*; as discussed 
in section 14. 1[ the contribution of this additional term can be neglected in the planar approximation 
of this study. 

The disturbing potential (i.e. the potential to be added to the potential of a sphere) due to ( , on 
an external point P*(r*, 0* , ip*), is 

SUo = 2k ^f eo (3 cos 2 W Q - 1) = - kfG ™^ £ ° ( sin 2 0* (3 cos(2ci - a Q ) + l) - 2 cos 2 p) (24) 

where G is the gravitation constant, a — Co/2 is the angle between the major axis of the spheroid and 
the meridian passing by the point P* and 

k f = ^4 (25) 

is the fluid Love number. We prefer to use kf instead of A, moment of inertia of m with respect to its 
major axis, because kf is an adimensional quantity. We remind that for a homogeneous sphere kj = 1.5 
and that this value decreases as the mass of the body is more concentrated in its central parfl From 

8 In standard theories using Love numbers, often kf is substituted by some smaller value to take into account 
the rigidity of the body. Here, a substitution of this kind is not allowed since kf is just being used as a convenient 
substitute for A. 
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now on, since the prolatcness appears multiplying all terms, we may use the mean radius R instead of 
R e . The differences due to this change are of higher order and may be neglected. However, when dealing 
with the first term of eqn. (f2"3|). we have to remember that, in the prolate spheroid, R e ~ R(l + eo/6). 

To obtain the force acting on one mass located in P*, we have to take the negative gradient of Uo 
and multiply the result by the mass of the point. It is worth emphasizing that (r* ,6*, ip*) are spherical 
coordinates of one point external to the body while in the preceding equations (£, 9* 1 (p*) represented 
the spherical coordinates of one point on the surface of the body. 

Since we are interested in the force acting on M due to the tidal deformation of m, once the gradient 

def 

is calculated, we can substitute (r*,9*,ip*) by the coordinates of M: (r, ^,<p = vo + v). We remind 
that this identification cannot be done before the gradient is computed because, in SU , we have both, 
the coordinates (r,<p) of the body M, and the coordinates (r*,ip*) of the generic point P* (where the 
gradient is taken). We thus obtain, 

k f GMmR 2 e / . \ 
Fi = - — 4 (3 cos(2i; - 21 - a ) + 1 J 

F 2 = '" (26) 

2k f GMmR 2 e . . 
F 3 = t —— sin (2v -21- flo - 
or 4 

The component F% (force along the meridian) is zero because we are considering here only the planar 
case. The corresponding torque is 

Mi = 

2k f GMmR 2 e . . , 
M 2 = 3 -sm{2v-2i-(^) (27) 

M 3 = 0. 

For the other terms of C(<?*) we have similar expressions having just to pay attention that, in these 
terms, k£ appears added to the arguments. We can proceed in the same way as above because the 
operations done involve only geometric quantities. We thus have 

5U= - kfGm ^ 6k ( sin 2 6* (3 cos(25 + U - a k ) + l) - 2 cos 2 fl* ) (28) 

k=—N r 

where the Ck are the prolatenesses defined by Eqn. (|2"Tj) . The corresponding force and torque are 
Fl = ^3 -(3cos(2v + (k-2)e-a k ) + lj 

k=-N "' 

A 2k f GMmR 2 e k ./„,,„,„ N 
F 3 = —^a -sm(2v + {k-2)£-a k ) (29) 

k=-N T 

A 2k f GMmR 2 e k , „ x<1 N 
M 2 = V 1 3 sin (2v + (k-2)£-a k ). 

of 

k=-N 



4.1 The axial terms 

At last, we have to consider the term R", not yet considered, and the terms ^i? e efeCOS 2 9* subtracted 
from the parts of £ to complete the equation of the spheroids. Putting then together, we obtain 



^axC = \R e ( 6k ~ e ~ 7 ' / £ P e7 **) cos 2 9*. (30) 



k=-N 



r = R + (5 ax C is the equation of a spheroid with symmetry axis perpendicular to the orbit and whose 
prolateness is given by the bracket in the above equation. As is well-known, the resulting field is axial. 
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The force on M is central and its contribution to the torque is null. It will be important in the general 
non-planar problem because it will contribute for the precession of the axis of the body. Because of 
its sign, it will counteract the effects due to the oblateness of the body (not considered in the present 
study). It will contribute short-period variations in the semi- major axis and eccentricity, which will be 
averaged to zero over one orbit. Hence, they can be neglected in the planar case. 



5 Rotation of close-in companions 

We use the equation Cfi = M2 (= — M z ) (see FRH Sec. 7-8). The time average of Si over one period 
is 



< Si > = ^ / M 2 dl 
2-kC J q 

Hence, using the approximation A ~ C in kj and simplifying: 



• A5GM 2 R 3 
< Si > = 



16ma' 



f 



1 2 1 4 13 6 \ . „ /49 2 861 4 21975 6 \ . „ 

/289 4 1955 fi \ . „ 1 fi . „ 714025 fi . „ 

e 4 e b sm2cr_2 H e b sm2cr 3 H e b sm2a-_ 3 

V 4 6 J 2304 J 2304 3 

It is worth mentioning that only squares will contribute to the average and the above result may 
be written as 

• A5GM 2 R 3 A . , . „ 
<f2> = 16^^ E El k (e)sm2a k . (32) 

k=-N 

In order to have an explicit equation in terms of the relaxation parameter 7 and the involved 
frequencies, the definitions given by eqs. (|19p may be introduced into the above equation through the 
trigonometric relation s'm2X = 2tanX/(l + tan 2 X), that is 

. „ 2^f(v + kn) ,„„. 

7 2 + (^ + fcn) 2 v ; 

The resulting expression can be used to study the tidal despining of close-in companions and/or 
central bodies. 



6 Synchronization. Spin-orbit resonance 

The immediate consequence of eq. pip is that the synchronous rotation is not a stationary solution of 
the system when the orbital eccentricity is not zero. Indeed, introducing eq. (|33[) and making v = 0, 
there results, in the first approximation, 



< n > 



135GM 2 R 3 nje 2 

~ ■ — . (34) 

u=0 2ma 6 (n 2 +7 2 ) 



The equality to zero is not possible if je ^ 0. In the synchronous state, the torque is positive, meaning 
that the rotation is being accelerated by the tidal torque. The stationary solution can only be reached 
at a supersynchronous rotation. Indeed, solving the equation < fl >= 0, we obtain 

6n 7 2 2^0 2 226n 6 + 1453n 4 7 2 + 28n 2 7 4 + 7 6 4 , n( 6 , 

S2 = n^ -e +3n 7 — — — e + C(e u ). (35) 

n 2 +7 2 ' 8(n 2 +7 2 ) 3 (4n 2 +7 2 ) y ' y ' 

The result corresponds to a supersynchronous rotation. However, at variance with the standard 
theories, the stationary rotation speed is not independent of the body rheology. It depends on the 
viscosity 77 through the relaxation factor 7. 
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In the quasi-inviscid limit, T] — > 0, then 7 ^> n and 7 2+ n 2 — 1- We then obtain 

Aim ^ n(l + 6e 2 + ^e 4 ). (36) 
o 

Thus, in the quasi-inviscid limit, the result is the same obtained with Darwin's theory when we neglect 
the differences in the response factors ki and assume that the ad hoc lags of tide components with 
equal frequencies are equal (see Laskar et al. 2004; FRH, sec. 9) 



7 Energy dissipation 

The rate of the work done by the tidal forces is W or b = Fv. 

The components of v in the adopted 3D spherical coordinates are 

nae sin v 

Vl 



r 

The result, time-averaged over one period, is: 



v 2 = (37) 

na 2 \/l — e 2 
v 3 = - 



• 3kfGM 2 R 5 i 
< W > orb = ' 



4a 6 



(1 - 5e* + f - § *) sin2. + (^ - ^ + ^ sin2 , 1 + 



(38) 

,147 2 2583 4 65925 fix ,289 4 1955 Bs e 6 3570125e 6 

( — 6 "^2- e +^12- e)sm2(7 - 1 + ( — 6 " ^ e } Sm2a - 2 - 4^08 Sm 2(73 + ^60^ Sm2(7 - 3 . 

or 

3fc fGM 2 R?n x — ■* /„ , \ „9 . . c„_\ 
< W >orb= j 6 Y ( 2 ~ ^./.(e) sin 2(7*. (39) 

a fc=-iV 

In the pseudo-synchronous stationary rotation, we may use v as given by eq. (|35p . Hence, using the 
sin2cr/ c values given in eq. (1331) and neglecting terms of order higher than 0(e 2 ), 

75k f GM 2 R 5 ne 2 -in , 

< W > orb (stat) ^ L— __. (40) 

In addition, we have to consider the work done by the tidal torque on the rotating body: < W > ro t= 
CO < SI >, that is, 

• 3k f GM 2 R 5 f2 A ~ . , . „ 
< W > rot = f -^- 6 Yl El k (e)sin2a k . (41) 

k=-N 

(Since W oc SI, the work associated with the rotation of the body vanishes when it reaches the 
stationary state.) 

The rate of the mechanical energy released inside the body is 

<E>=-(<W> olh + <W> Iot ) >0. (42) 

GmM 2a 2 

From eq. ((38l) . since W or b = ~ , we obtain a = — — — Worbi i.e. the secular variation of the 

2a GmM 

semi-major axis 

3kfMR 5 n J^-y ,„ ¥ ,_o , , . „ , 

< a >= 1 Y, ( 2 " k)El k (e)sm2a k . (43) 

k=-N 

The interpretation, neglecting the terms in e 2 , is easy. If the body is rotating faster than the orbital 
motion, then v > 0, cto > 0, and d > 0. The tide in m causes the bodies M and m to recede one from 
another. Otherwise they are falling one on another. 
Two approximations of this formula are useful: 
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(i) The free rotating approximation 



and 



3kfMR 5 n~/is 

< a >ft ™" 17 2^ 2^ ( 44 ) 



(ii) The pseudo-synchronous approximation 



75k f MR 5 ne 2 jn 

<a> ^ t^r^- (45) 



7.1 The quality factor of standard theories 



The quality factor Q is a parameter originally introduced to characterize damped oscillators. It ex- 
presses the quality of the oscillator in keeping free oscillations alive. It is proportional to the proper 
frequency of the oscillator and vanishes when no elastic force is acting. Its extension to forced oscil- 
lations is not done without ambiguities and we do not use it in this theory. Nevertheless, the quality 
factor Q is widely used and, in the applications, we need to know how to express it in terms of the rheo- 
physical parameters used here. However, it is worth emphasizing that the formulas given in this section 
are obtained by mere comparison of some equations of this theory with their equivalents in standard 
Darwin theories and are not valid out of the particular conditions in which they were established. 

We recall that, in standard theories, the quality factor Q and the tidal Love number k 2 cannot be 
separated one from another; in the following equations, k 2 is the tidal Love number and kf is the fluid 
Love number. 

In standard theories, two different definitions of Q are used, one when the body is free rotating and 
another when it is trapped in a stationary pseudo-synchronous state. 



7.1.1 Bodies in free rotation 



In classical theory, in this case, we have 



< W >c 



3k 2 GM 2 R 5 n 
2a 6 Q 



(see FRH eq. 48) where Q is the inverse of the lag of the semi-diurnal tide (so). Comparing to the 
eccentricity- independent term of eq. (|38[) . we obtain the equivalence formula (valid only for small 
eccentricities) : 



Q = 



k 2 (7 2 + 



1 



i -i 



■ sin 2<7i 



X + 



(46) 



where we have introduced x — ~ (x is the frequency of the semi-diurnal tide in units of 7). 

It is important to note that Q goes to infinity when \ (or v) goes to zero. This is so also in the 
standard theories and is just a consequence of the inadequacy of the quality factor Q to measure 
dissipation. One may note that the dissipation itself, given by eq. (|39l) . is not singular for v = 0. 



7.1.2 Bodies in pseudo-synchronous stationary rotation 
In standard theories we have, in this case, 

75k 2 GM 2 R 5 ne 2 9k 2 GM 2 R 5 ne 2 
<W> ^ 8^Q 8^ 65 

(see FRH eq. 51) where, now, Q is the inverse of the lag of the monthly /annual tide (£2)- In FRH, we 
have considered the lag of the radial tide (£5) as equal to the lag of the monthly /annual tide (£2) since 
both have the same period in the case of a stationary rotating body. If we compare only the first part 
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Fig. 4 The low-eccentricity equivalent of the quality factors (in units of fe/fc/) as functions of the frequency 
(in units of 7). For rotating bodies, \ — v H (solid line); for bodies in stationary rotation, \ = 71/7 (dashed 
line) . 



of the above equation to eq. (|40p . we obtain again eq. f4"6")l . However, when the complete equation is 
considered, there results 



k 2 (7 2 + n 2 ) 



k 



7 



yn 



28 k 2 
25 ~k~ s 



— sin 2<7i 
2 



28 fc 2 / 1 
25^7 V + X 



(47) 



where now x = ^ (x is t ne frequency of the monthly/annual tide in units of 7) 

This duality in the actually used definitions of Q in standard theories is a big nuisance. In the case 
of planetary satellites, eccentricities are low and either the system is rotating or nearly synchronous 
and we may consider the two cases separately. However, in the case of exoplanets, eccentricities are 
often high and the choice of one of the two formulas to determine Q is a problem. Indeed, in the 
standard approach, if the eccentricity is large, the stationary rotation may have a period much smaller 
than the orbital period and the dissipation due to the semi-diurnal tide will not vanish as in true 
synchronous companions. As a consequence, tidal components with frequencies v and n will contribute 
to the energy dissipation on the same foot making impossible to privilege one of them to define one 
quality factor. 



8 Circularization 

The variation of the remaining elements can be obtained straightforwardly using Gauss equations (see 
Bcutler 2005, Sec. 6.3.5). As discussed in FRH (section 18.1), in order to take into account correctly 
the reaction on M of its tidal action on m, the accelerations R' , S' , W of those equations need to be 
multiplied by (M + m)/m or, equivalently, by n 2 a 3 /Gm. With the forces calculated in Section |H we 
thus get 



3k f MR 5 ne 

< e >= o r 

oma° 



, 21 , . 3299 BN . „ ,1 1 2 35 4 175 fiN . „ 
(1 - "J + ^ - -576" e ') Sm2 ^ + (4 - I6 G - 768 £ " 18432 6 ) 



(48) 



, 49 1253 2 50311 4 508651 fiN . „ , 289 2 10421 4 310463 fiN . 

( 1 e 2 e 4 H e 6 ) sin2cr_i + ( e 2 H e 4 e 6 ) sin2<T_ 2 + 

V 4 16 256 2048 ' V 2 12 144 ; 

/ 1 4 17 as . „ ,714025 4 105299675 BN . „ e 6 . „ 284089e 6 . „ 

e -| e sm 2(73 — e e sin 2<t_t H sin 2(74 sin 2er_4 

v 768 18432 ; J v 768 18432 ; 6 144 64 



or 

ZkfMR? X 

< e >= - 



8ma 5 e 

k=—N 

Taking into account eqs. (1331) . we may also write 



(2-v/l - e 2 - (2- fc)(l - e 2 ))j5| ifc (e) sin2cr fc . (49) 



3kfMR 5 ney ( . v (v — n) (v + n) , , . 

16ma 4 V 7 2 + v 2 j 2 + (v - n) 2 j 2 + (v + n) 2 ' 
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In the case of a stationary or near-stationary rotation, v — 0(e 2 ) and the above equation is reduced 

to 

75 k f MR 5 ne-f n , . 

<z>--^~ 5 1 2 ^ 2 - (51) 

8 mo? 7^ + n z 

9 Dissipation parameters in stars, planets and satellites 

In this section we determine the values of 7 for several Solar System and extrasolar bodies. We use for 
that sake the values published in the literature usually obtained using standard tidal evolution theories. 
One problem common to most of the given examples is that the inversion of the equivalence formulas 
relating 7 to Q has two solutions. The choice of one of the two solutions is done after comparing the 
values of the equivalent uniform viscosity in each solution. 



9.1 Io 

The tidal evolutions of the Galilean satellites of Jupiter are among the best studied in our Solar 
System. From the satellites accelerations, Lainey et al. (2009) have determined the dissipations of Io 
and Jupiter. For Io, they have found k 2 /Q = 0.015 ±0.003. Introducing this value in the formulas given 
in section l7Tl we obtain 7 = 4.9 ± 1.0 x 10~ 7 Hz. It is worth mentioning that this result is independent 
of the individual values of k 2 and Q, which are not well known. The calculated value of 7 depends only 
on the value of k 2 /Q and on the moment of inertia (0.378 mR 2 ). 

With the results given in Sec. 151 (eq. I3"5l) . we obtain for the synodic rotation period (a.k.a. length of 
the day) P syn = 3300^Jqqq yr. We may compare this value with the minimum value 1400 yr determined 
by Milazzo et al. (2001) from the comparison of Galileo and Voyager images taken 17 yrs apart. 

The equivalent viscosity corresponding to this determination may be obtained using eq. The 
result, 1.2 ± 0.3 x 10 16 Pa s, is in good agreement with the value used by Segatz et al. (1988): 2 x 10 16 
Pa s, in their models for Io's tidal dissipation. 

We note that the adopted value of 7 corresponds to \ — 94, that is, to the ascending branch of the 
curve shown in fig. 2] (Darwin's theory corresponds to the descending branch). Therefore, we are in the 
regime proposed by Efroimsky and Lainey (2007) in which Q increases with the frequency of the tide 
component. The solution corresponding to Darwin's regime gives 7 = 0.003 and t] = 1.7 x 10 12 , which 
is 2 — 3 orders of magnitude below the known viscosity of ice and silicates (resp 1.5 x 10 14 Pa s and 
3 x 10 15 Pa s cf. Sotin et al. 2004). 



9.2 Europa 

The basic information we have on Europa is its non-synchronous rotation. The synodic period is 
constrained by the lower limit 12,000 yr obtained from the comparison of Galileo and Voyager images, 
and the upper limit 250,000 yr, obtained by comparing the present position of some cycloidal cracks 
with the longitudes at which their shapes should have been formed (see Greenberg et al. 2002). Using 
the results given in section [J51 we obtain 7 = 1.8 - 8.0 x 10~ 7 Hz and k 2 /Q between 0.01 and 0.045 (if 
we adopt k 2 = 0.26, we obtain Q in the range 6-26). The equivalent uniform viscosity corresponding 
to these results is 77 = 4 — 18 x 10 15 Pa s. 

The same indetermination discussed above occurs here since the excess of rotation speed is approx. 
proportional to x 2 + x -2 - As before, we have to look to the viscosity to decide between the two 
mathematically possible solutions. The solution corresponding to Darwin's regime gives 7 = 0.0005 — 
0.002 and r\ = 1.3 — 6 x 10 12 Pa.s, which is less than expected, but not so expressively as in other cases. 

We note that the indetermination is as more difficult to solve as Q is small. 



9.3 The Moon 

The quality factor of the Moon has been determined for two tidal frequencies. For the monthly tide, 
Q = 30 ± 4 and for the annual tide Q ~ 35 (Williams et al., 2005; Williams and Boggs, 2008). The 
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values of Q for the monthly tide (combined with the very low ki of the Moon, 0.0301 cf. Williams 
et al.) gives 7 = 2.0 ± 0.3 x 10 -9 Hz. The synodic period of the Moon corresponding to the above 
determinations of 7 is larger than 5 Myr and indistinguishable from a true synchronous spin-orbit 
resonance. 

The problem with this determination is the value found for the uniform viscosity: rj = 2.3±0.3 x 10 18 
Pa s, 4 — 5 orders of magnitude smaller than the values reported in the literature, which refer to the 
solid lithosphere. This seems to be one more indication in favor of the role played by a plastic lunar 
asthenosphcrc in tidal dissipation and is in agreement with the dramatic decrease of the seismic Q 
below the deep moonquake source region indicating the presence of a partial melt below this depth 
and that, likely, most of the solid body dissipation in the Moon occurs below 1150 km (see Wieczorek 
et al. 2006). The alternative solution, corresponding to Darwin's regime, gives viscosity values yet 
smaller. The adoption of a more realistic value for the viscosity leads to a much smaller 7 and a much 
larger synodic period, even when a permanent equatorial asymmetry of the Moon is neglected. The 
resulting 7 is small enough to compensate the fact that n and e may have been larger in the past and 
to give a result consistent with the rotation indicated by the distribution of the craters on the Moon: 
a synchronous attitude lasting since the formation of the last great basin, 3.8 Gyrs ago (see Wieczorek 
and Le Feuvre, 2009). However, a larger viscosity would also imply a larger Q, in contradiction with 
the above-mentioned determinations^ 

Let us add that the annual tide cannot be studied with the Keplerian model used in this paper. It 
is clearly related to the motion of the Earth-Moon system around the Sun and a perturbed model is 
necessary to interpret it. 

9.4 Titan 

Cassini radar observations of Titan over several years show that the present-day rotation period of 
Titan is different from synchronous and correspond to a shift of ~ 0.12° per year in apparent longitude 
(Stiles et al. 2008, 2010). This result, if only due to tidal torques, would mean a very large dissipation 
(Q < 10). Values of Q are not reported in recent literature, however, from internal structure studies, 
Tobie et al. (2005) conclude that because of convection in the outer ice layer, we have now Q ~ 50 (but 
Q > 300 during almost the whole satellite past evolution). Adopting this value and ki in the interval 
0.32-0.39 (cf Sohl et al. 2003), we obtain 7 = 2.9 ± 0.2 x 10~ 8 Hz, 77 = 1.1 =b 0.1 x 10 17 Pa s and a 
synodic period P ~ 200,000 years. As the Moon, Titan behaves much like a solid body. This favors 
the interpretation of the measured shift of Titan's crust as due to seasonal effects (e.g. the interaction 
between the crust and the atmosphere; see Tokano and Neubauer 2005; van Hoolst et al., 2008). 

Let us also mention that a re-analysis of the Cassini data by Meriggiola and less (Meriggiola, 
2012) has not showed discrepancy from a synchronous motion larger than 0.02° per year. This more 
conservative result does not disagree with those of the analysis done above. It only sets less strict limits 
as 7 < 1 x 10~ 7 Hz and Q > 15. 

9.5 Jupiter 

From the acceleration of the Galilean satellites, mainly Io, the dissipation in Jupiter has been deter- 
mined to be given by k^jQ = 1.102 ± 0.203 x 10~ 5 (Lainey et al. 2009). This value corresponds to 
X = 1.1 ± 0.2 x 10~ 5 and 7 = 23 ± 4 Hz. The equivalent uniform viscosity is r\ = 4.7 ± 0.9 x 10 10 Pa s. 
Jupiter's tides are in Darwin's regime (the other alternative would need a viscosity as high as 10 20 Pa 

s) 

9.6 Saturn 

From the maximum possible past evolution of Mimas, Meyer and Wisdom (2007) obtained for Saturn, 
the limit Q > 18, 000. If we use Saturn's Love number, ki = 0.341, we obtain x < 2.4 x 10~ 5 (Dar wm s 
regime) and 7 > 13 Hz. The equivalent uniform viscosity is r\ < 1.5 x 10 10 Pa s. 



9 In the Efroimsky-Layney regime, Q increases with the relative frequency \ of the tide component and then 
increases when the critical frequency 7 decreases. 
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Table 1 Summary of the values adopted and/or obtained in this section. 



Body 


7 (Hz) 




2tt/7 


r\ (Pa s) 


Q equivalent 


Moon 


2.0 ±0.3 x 10" 


9 


36,000 d 


2.3 ±0.3 x 10 18 


30 ±4 


Titan 


2.9 ±0.2 x 10" 


8 


2500 d 


1.1 ±0.1 x 10 17 


~ 50 


Solid Earth 


0.9-3.6 x 10" 


7 


200-800 d 


4.5 - 18 x 10 17 


200-800 


Io 


4.9 ± 1.0 x 10" 


7 


730 d 


1.2 ±0.3 x 10 16 


50 - 80 


Europa 


1.8-8.0 x 10" 


■7 


90-400 d 


4 - 18 x 10 15 


6-26 


Neptune 


2.7-19 




< 2 s 


1.2 - 4.8 x 10 10 


9,000-37,000 


Saturn 


> 7.2 




< 0.9 s 


< 15 x 10 1(1 


> 18, 000 


Jupiter 


23 ±4 




~ 0.3 s 


4.7 ±0.9 x 10 10 


~ 36,000 


hot Jupiters 


8-50 




0.1-0.8 s 


5 x 10 10 - 10 12 


2 x 10 5 - 2 x 10 6 


solar-type stars 


> 30 




< 0.2 s 


< 2 x 10 12 


> 2 x 10 6 



9.7 Neptune 

Founded on previous studies of the Neptunian satellites, Hamilton (2009) proposes the value Q/kz = 
4.5 x 10 4 (with an error factor 2). If we adopt the Love number, hi — 0.41 (Durda, 1992), there follows 
7 = 9.4 Hz and r\ — 2.4 x 10 10 Pa s. The factor 2 of incertitude in Q is reproduced in the values of 7 
and r\. 



9.8 Solid Earth 

The tidal dissipation in the solid Earth was estimated from satellite tracking and altimetry (Ray et 
al. 1996) as Q = 370 with an error factor ~ 2 (lag angle 0.16 ± 0.09 degrees). Considering the Love 
number k2 = 0.46, there follows \ — 770, 7 = 1.8 x 10~ 7 Hz and rj — 9 x 10 17 Pa s. All these values 
have error factors ~ 2. The viscosity found is several orders of magnitude smaller than the viscosity of 
the Earth's lower mantle (see Karato, 2008). This uniform value does not make a distinction between 
the various parts of the solid Earth and is expected to be smaller than the viscosity of the mantle. 
The result, however, allows us to solve the two-solution indetermination in the inversion of Q(x)- For 
instance, in the present case, the root corresponding to a Darwin's tide regime (7 ~ 0.1 Hz) leads to 
viscosity values 6 orders of magnitude smaller than the above given one and can be excluded. The low 
viscosity found is an indication that the given 7 is too large and that we should look for a creeping 
law leading to more dissipation in high frequencies, as it happens when the Andrade model is plugged 
in the modified standard theory as proposed by Efroimsky (2012). 



9.9 Hot Jupiters 

A first assessment of the values of 7 for hot Jupiter may be done using the values found by Hansen 
(2010) from an analysis of the survival of some short-period exoplanets of mass ~ 0.5Mj up . The 
comparison of our results to the formulas used by Hansen adapted to the case of a planet trapped in 
a stationary rotation (i.e. pseudo-synchronous), gives 

^=421^ (52) 

where a p is the planetary dissipation parameter used by Hansen. 

Using his mean results for WASP-17 b, Corot-5 b and Kepler-6 b, transiting planets whose radii 
have been determined, we obtain values 7 in the range 8-50 Hz, the smaller value corresponding to 
the bloated WASP-17 b (radius ~ 2i?j up ). The results are also sensitive to the moment of inertia of 
the planets (via kf) and were obtained using A > O.lmR 2 . If the central concentration is yet larger 
(i.e. if A is smaller), 7 may be smaller. 

Another possibility is to use the known value of Jupiter's Q and some scaling laws : (i) If n <C 7, 
Q scales with the period of the main tidal component (see section ITTTj) ; (ii) Q scales with R~ 5 (see 
Eggleton et al. 1998; Ogilvie and Lin, 2004, Hansen, 2010). For instance, we may consider one hot 
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Fig. 5 Simulation of the tidal evolution of the orbit of planet CoRoT 5b using 7 = 200 Hz (black) for 
initial eccentricities 0.05, 0.09 and 0.18. The corresponding results with the Darwin-Mignard approach using 
Q — 3.4 x 10 6 are also shown (blue). Vertical lines show t = and the age of the star range (5.5-8.3 Gyr). 
Horizontal lines: Observed eccentricity at t = (limits of the error bar c/. Rauer et al. 2009). 



Jupiter with a mass of 2-3 Mj up in a 5-day orbit. The period of the main tide raised on it is 29.4 times 
the period of the semi-diurnal tide of Jupiter (raised by Io) and planets of this mass range have radii 
1.2 ± 0.2i? Jup . With these data we obtain Q ~ 420, 000 and, then, 7 = 15 Hz. 

The very existence of hot Jupiters around old stars in significantly non-circular orbits is an im- 
portant test for tidal theories. Indeed, all general mechanisms responsible for important eccentricity 
enhancement are related to events expected to occur in the early stages of the formation of the system 
(see Malmberg and Davies, 2009). Therefore, in older systems, if eccentricities were not damped to 
zero, the variations due to the tidal evolution may have been below some limits. We note that the 
existence of several planets in these conditions^ pl a y s against hypothesizing that exceptional sources 
of significant enhancement may have existed in the recent story of each one of them. However, we 
emphasize that the two parameters considered in this analysis, age and eccentricity, are of difficult 
determination. 

We have studied some of the C0R0T hot Jupiters in elliptic orbit. On one hand, transiting planets 
have better determined eccentricities^!, and, on the other hand, their ages were determined with some 
confidence. 

Fig. [5] shows the tidal evolution of the hot Jupiter CoRot-5 b calculated using the approach devel- 
oped in this paper (black lines) and the standard Darwin theory (blue lines), respectively. In order to 
avoid well-known truncation errors associated with expansions, we have used 7V=150 in the expansions, 
which allows dealing with eccentricities in the range of the solutions shown (e < 0.85) without serious 
truncation errors (see the Appendix). The evolution predicted by Darwin's theory was simulated by 
numerical integration of one 2-body model using Mignard's expression for the tidal force in closed 
form. 

Fig. [5] shows the good agreement of the rheophysical theory of this paper and Darwin's standard 
theory in the study of the evolution of giant planets. The solutions shown in figure [5] correspond to 
7 = 200 Hz and its equivalent Q = 3.4 x 10 6 , respectively. 



10 e.g. CoRot-5 b (M = 0.47A'/. Jup , e = 0.09±g;^), CoRoT-12 b (M = 0.9M Jup ,e = 0.07±g;^) and CoRoT-23 
b (M = 2.8Af.j up , e = 0.16 ± 0.07) 

11 The true longitude must be 90 degrees at the minimum light of the transit. This constraint added to 
the radial velocity measurements, allows a better determination of the eccentricity and the longitude of the 
pericenter than in the case of non-transiting planets where these parameters are to be determined based on 
hard-to-measure asymmetries of the radial velocity curve. 
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9.10 Hot super-Earths 

The few known transiting planets in the 1-10 Earth mass range are in circular orbits and the memory 
of their past evolution is erased. There is only one case among the currently known ones that may give 
us some information. It is 55 Cnc e for which e = 0.057^q'q^[. It belongs to a somewhat hierarchized 
system of 5 planets whose past evolution may be simulated from ad hoc initial conditions able to bring 
the system to its present situation. 

Simulations with the standard theory (Mignard's torque) starting with arbitrary excited eccentric- 
ities showed that dissipation in 55 Cnc e drives the innermost planets (55 Cnc e and 55 Cnc b) to a 
stationary solution with aligned pericenters (secular mode I of Michtchenko and Ferraz-Mello, 2001) 
and the eccentricity of 55 Cnc e falls to ~ 0.003 in a time much shorter than the age of the system 
(which is 10.2 ± 2.5 Gyr). If Q > 5500 (dissipation factor of the current "annual" tide), the damping 
to the equilibrium center of the secular dynamics is much slower allowing 55 Cnc e to be found at a 
larger eccentricity now. 

The only other possible guess comes from the viscosity of CoRoT-7 b estimated by Leger et al. 
(2011): i] > 10 18 Pa s. Using the physical data of that planet, we obtain 7 < 5 x 10 -7 Hz and Q — 100. 
This result means that for CoRoT-7 b, % > 1 and Q grows with the frequency (Efroimsky-Lainey 
regime) , one fact that should be taken into account in the simulations of the evolution of the system. 
If tt2 = 0.46 as for the Earth, then Q' = 3Q/2/c2 = 300, which is of the same order as the value 100 
adopted by Rodriguez et al. (2011) in the study of the tidal evolution of the CoRoT-7 planets. 



9.11 Solar type stars 

Dissipation values of solar type stars have been estimated by Hansen (2010) from an analysis of the 
survival of some short-period exoplanets. The comparison of the equations used in Hansen's evolution 
model to ours allows us to write: 

7 = 'f /Star (53) 

t3 - ft star cr star 

where er s tar is the stellar dissipation parameter used by Hansen. Using his mean result cr s t a r = 8.3 x 10 -64 
g^ 1 cm -2 s" 1 , we obtain for M and G stars, 7 ~ 3 — 25 Hz. This variation is mainly due to important 
dependence on the radius and the upper limit correspond to stars having half of the radius of the Sun. 
For one star equal to the Sun, the result is 7 = 7 Hz. 

An approximated calculation gives for the corresponding kinematic viscosity, 10 9 — 10 ms, 
which are much larger than the value 10 8 m 2 s used by Ogilvie and Lin (2007) in their models of 
tidal dissipation in stars. 

However, in the study of transiting hot Jupiters mentioned in Sec. 19.91 we have found that 7 < 30 
Hz often implies a too large exchange of angular momentum between the orbit and the star rotation. 
As a consequence, we needed to assume extremely low values for the star rotation in the past to be 
able to reproduce current observed values. 

The limit 7 > 30 Hz given in table 1 is more or less the same obtained by Jackson et al. (2011) from 
the analysis of the distribution of the putative remaining lifetime of hot Jupiters. We note that the 
corresponding limit for the viscosity is 30 times larger than that adopted by Ogilvie and Lin (2007). 



9.12 Binary stars 

We may use a sample of data on detached binary stars selected from those collected by Torres et al. 
(2010), to investigate their rotations. If we fit a law Ae 2 through the points shown in fig. [6l we obtain 
a coefficient close to 8. However, a further analysis shows that this value is strongly determined by 
the few points corresponding to e > 0.2; when these points are not included the coefficients falls to 
2. The concentration of the measured values around the synchronous rotation is clearly seen; we may 
also see that among the stars in elliptic orbits a majority shows rotation above the synchronous value 
(i.e. they are super-synchronous). However, the high dispersion of the points does not allow us to fix 
the value of A, nor to discard the coefficient corresponding to a standard solution (A = 6), which may 
prevail because of the low viscosity (and consequently high 7) of normal stars. 
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Fig. 6 Binary stars: Distribution of the angular velocity of rotation in function of the eccentricity. 
10 The elastic tide 

The shape of the body deformation due to the creep tide is given by eqn. (J2DJ) . After the transient phase 
(i.e. for jt 3> 1), only the forced terms matter and it is dominated by the semi-diurnal component 



the maximum of which is reached when 2cv — Co = 0, i.e. the angle between the vertex of this component 
to the sub-M' point is ctq/2. We remind that Co is a constant determined by the integration of the creep 
equation and not an ad hoc lag; it is not necessarily small. In addition, tanero is proportional to the 
semi-diurnal frequency v and inversely proportional to the relaxation factor 7 (i.e. it is proportional 
to the viscosity). 

In the case of an inviscid body, the response to the tidal action is instantaneous, 7 — > 00 and 
o"o — > 0; the tide highest point remains aligned with the mean direction of the tide raising body M (as 
shown in Fig. [J). 

However, when 7 <C i>, as in rocky planets and satellites, ero will approach 90° (in the rigid body 
limit, 7 = and ctq = 90°). This result is in contradiction with the observations. For instance, the 
observed geodetic lag of the Earth's body semi-diurnal tide is very small (0.16 ± 0.09 degrees cf. Ray 
et al. 1996). In addition, some authors (see Efroimsky, 2012) claim that in these bodies the actual lags 
do not obey the weak friction approximation where lag tangents are proportional to the frequencies, 
but have lags proportional to a negative power of the frequency. 

In order to conciliate the theory and the observed tidal bulges in the Earth, we have to assume 
that the actual tide is not restricted to the component due to the creeping of the body under the tidal 
action, but has also a pure elastic component. No matter how empirical this hypothesis seems to be, 
it explains well the observed behavior. 

Let this elastic component be defined at each point by its height over the sphere and be given by 
S( e i((j>* ,6*) — X(p((j>* ,9*) — R') where p is the radius vector of the equilibrium spheroid surface and 
A is a quantity related to the maximum height of the tide (see Sec ll0.2[) . For the Earth, for instance, 
A ~ 0.2, which is the ratio of the observed maximum height of the lunar tide (26 cm after Melchior, 
1983) to the maximum height of the equilibrium spheroidal figure (1.34m). 

The sum of the (local) heights of the elastic tide and of the main term of the creep tide is 



where, for the sake of simplicity, we have set i?2,o(e) = 1 and sin#* = 1 (equator). 

There is some similarity between this composite model and the model studied by Remus et al. 
(2012) and due to Zahn (1966). The elastic tide of this paper is, in principle, the same as Zahn's 




(54) 




(55) 
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Fig. 7 (a) Geodetic lag of the semi-diurnal tide as a function of <to. (b) Same as (a), but as a function of 
X = v/'y (x — tan (To), (c) Time evolution of the geodetic tide lag when the frequency of the semi-diurnal tide 
crosses and the tidal bulge changes of side with respect to the sub-M point, (d) Maximum height of the tide 
in units ^Re p . The blue lines correspond to A = 0.2 (Earth). 



adiabatic tide. However, the creeping tide is very different from Zaun's dissipative tide. The physical 
setting of the two models is not the same and the results are different; for instance, the creeping tide 
is not in quadrature with the exciting potential, as Zahn's dissipative tide. 



10.1 The geodetic lag 

The maximum tide height (i.e., the maximum of 5£) is, now, no longer reached at a — <7q/2 as the 
creep tide, but at 

a = -s (56) 

where 

sin2ero 

e = arctan- — — . (57) 

l + 2A + cos2o- 

This function is shown in fig. Efa). We see that, as far as A 7^ 0, Eq —> when ctq — >• f, that is, 
when 7 — > 0. As a bonus, we also have near the rigid limit (i.e. near Co = ■§) £ o decreasing when uq 
increases, that is, when the frequency v increases. This is exactly the behavior that is being advocated 
by Efroimsky and collaborators (Efroimsky and Lainey, 2007; Efroimsky and Williams, 2009; Castillo- 
Rogez et al. 2011; Efroimsky, 2012) for the Earth and the planetary satellites. This is best seen if we 
plot these curves using as independent variable the semi-diurnal frequency \ (in units of 7). See fig. 
EJb). We see that for \ small, £0 grows with % (as in Darwin's theory). However, in this model, it 
only grows up to reach a maximum and then decreases. We may compare these curves to the curve 
Q(x) presented in fig. 0] and remind of the popular use of the relation Q — l/e in Darwinian theories. 
However, the matching is very imperfect: The minimum of Q happens for x = 1 while the maximum 
of £q happens for values of x larger than 1, which increases indefinitely as A tends to zero. 
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It is worth recalling that one of the difficulties created by the assumption that the actual tide lag 
is proportional to a negative power of the frequency happens when the frequency changes of sign. If 
no additional assumption is done, we just have a singularity with the tide lag tending to infinity or, 
at least, abrupt jumps between positive and negative values. This is not the case with the solution 
that results from the superposition of the elastic and creep tides. In this case, the transition from 
one side to another is smooth. The tide angle increases when the frequency decreases up to reach a 
maximum; after that point it quickly decreases up to cross zero with a finite derivative. The behavior 
in the negative side is just symmetrical (see fi£0c) ). 

In addition, the agreement of fig Etc) with the red curve of Figure I of Efroimsky (2012) is note- 
worthy. Efroimsky's curve corresponds to the left half of fijQc). The only differences are the reversed 
direction (fig. 0(c) is drawn assuming that the frequency is decreasing as the time increases) and the use, 
by Efroimsky, of logarithmic scales. The two figures show exactly the same behavior notwithstanding 
the fact that they arise from two completely different models. 



10.2 Tide maximum height 

The maximum height of the creep tide (after the transient phase is over) can be determined as a 
function of the semi-diurnal frequency and the relaxation factor. It is given by the value of h = ( — R' 
at a = a max , that is 

h m ax = C(a m ax) - R' ■ (58) 



Hence, for 3> 1 and 9* — 90° (i.e. at the equator of m), in this approximation: 

1 T;Re P iRjep 

h max = -Repcosao = — f „ = , = (59) 

2 y/l + tan 2 ct Vl + ' 



This result is the same obtained by Darwin (1879) for the tides of a viscous spheroid, when inertia is 
neglected. 

In the case of one solid body, 7 = and the creep tide height vanishes. In the case of an inviscid 

body, 7 — ► co, uo —> 0, and we get the value h max = ^Rep- 
It is worth noting that the latest results appear clearly in the solution of the differential equation 

(eq. [7]) when it is written as 

C = Ce~ 7t +R'+ ^Re' p cosct cos(2a - a ). (60) 

When an instantaneous elastic tide is added to the creep tide, the maximum height of the geodetic 
(or composite) tide is the value of the function <5£ eq (eqn. I55|) at its maximum: 

1 

2' 



Lax - -i?e pV / A 2 + (l + 2A) cos 2 a . (61) 



Hence, it is almost unchanged when the viscosity is low (v <C 7 i.e. Co "C 1), but the difference becomes 
significant in the case of viscous bodies, when 00 is high. The relative value of the maximum height of 
the actual tide is shown in fig. E[d). In that figure, the unit is the maximum height of the equilibrium 
spheroid (^Re p ). One may note that when 7 <g; z^, 00 — >• 7r/2, the height of the creep tide tends to zero 
and the maximum height of the composite tide is the maximum height of the instantaneous elastic 
tide: ^XRe p . 

It is important to emphasize that the frequency-dependent height of the tide has not been taken 
into account in the majority of modern tide theories (Jeffreys, 1961; MacDonald, 1964; Kaula. 1964; 
Singer, 1968; Mignard, 1979; Hut, 1981; Murray and Dermott, 1999; Laskar and Correia, 2004). In these 
theories, Love's elasticity theory is used to calculate the potential of the tidally deformed body, which 
is proportional to the frequency- independent Love number &2 (or, in higher-orders, made of spherical 
harmonics proportional to the frequency- independent kj, j — 2,...). They consider the tides on an 
elastic body and, to take into account the imperfect elasticity, just delay the potential introducing by 
hand a phase lag. 
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10.3 Dynamical consequences of the elastic tide 

This empirical introduction of an elastic tide is necessary to have the theory conform to the measure- 
ments of Earth's bodily tides. It is however important to stress the fact that, being elastic, it does not 
present a lag. The major axis of the prolate spheroid corresponding to the elastic tide is permanently 
oriented towards M . 
It then follows, 

4k f GMmR 2 \e' 

F lcl = i — 4 £ = F 3el = 0. (62) 

The force due to the empirical elastic tide is radial and its torque is equal to zero. Therefore it does 
not contribute to the variation in the rotation discussed in section [5] Its contribution to the dissipation 
also vanishes since the term added to Fv is proportional to r sinn whose time average is zero. So, 
its introduction does not affect any of the results presented in previous sections 



11 Conclusions 

We will be brief in these conclusions. The rheophysical theory presented in this paper is yet incipient 
to justify a long discussion. However, even at this early stage, some results are noteworthy. 

The first result concerns the problem that served as motivation for this investigation: the so- 
called pseudo-synchronous stationary rotation. In this theory, it is given by a law that depends on 
the viscosity of the body. In the limit 7 — > 00 (or rj — > 0), the body behaves like a perfect fluid and 
the pseudo-synchronous stationary rotation happens exactly as in classical Darwin's theory. However, 
when 7 <C n (e.g. in rocky bodies), no matter if the eccentricity is large or small, the stationary rotation 
is close to the true synchronous rotation. In the case of rocky bodies, the monthly/annual tide (due to 
eccentricity) does not create the strong torque responsible for the super-synchronicity of fluid bodies. 

In the case of natural satellites, the reophysical theory allows us to obtain the excess of the ro- 
tation period and to determine the "length-of-day" (a.k.a synodic rotation period). The results are 
in agreement with the observed values. They were obtained without adding any torque due to some 
"permanent" equatorial asymmetry, which may indeed exist, but is not necessary to counterbalance 
the tidal torque. 

v 7 _i 

The second result concerns dissipation. Energy dissipation appears proportional to ( — I — ) . The 

7 v 

maximum dissipation of a tide component is reached when its frequency equals 7 and decreases sym- 
metrically when v is different of 7 no matter if larger or smaller. When v <C 7 the energy dissipation is 
proportional to v. This is what happens in giant planets and stars. When v 3> 7, the energy dissipation 
is inversely proportional to v. This is the behavior of dissipation in rocky bodies as extensively dis- 
cussed by Efroimsky and Lainey (2007), Efroimsky and Williams (2009), Castillo- Rogez et al. (2011) 
and Efroimsky (2012). We have taken some care writing these sentences to avoid privileging 7 or v. 
In fact, the relaxation factor 7 plays a role of critical frequency and one body may behave in different 
ways under the action of two tidal components of different frequencies if one is larger and the other 
is smaller than the critical frequency 7. In table 1 we have included in one column the inverse of 7, 
to give a clearer idea of the location of this bifurcation, which, in some cases, is not very far from the 
periods of the actual stresses acting on the body. 

At this point, I would like to stress that we have not used the quality factor Q in the present theory. 
In standard theories, the quality factor Q is an ambiguously defined parameter. It is defined using the 
semi-diurnal tide in the case of a freely rotating body, but using the monthly/annual tide in the case 
of a pseudo-synchronous companion. In classical applications, these two cases are very distinct and 
we can handle the two different definitions. However, in the case of high-eccentricity exoplanets, this 
separation no longer occurs. The dissipation is shared in comparable parts by the semi-diurnal and 
the monthly/annual tide and we get different values of Q following we consider one tide component 
or another. The low-eccentricity equivalence formulas relating Q to 7 and the values listed in table 
1 only appear in this paper because it is important to have a bridge between standard theories and 
this new one. However, it is necessary to stress the fact that these formulas are just numerical bridges 
valid for low eccentricities and for the rotational states assumed to establish them. Strictly speaking, 
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a universal relation between Q and 7 does not exist. We also note that the tidal Love number ki has 
not been used in the equations of the creeping tide. 

Section [9] presents a short inventory of bodies in the Solar System and extrasolar. For each of them, 
the values of the relaxation factor 7 and the uniform equivalent viscosity are derived on the grounds of 
the results obtained for them with standard theories. We stress the fact that the results shown allow 
the studied bodies to be divided into two groups. One group, including the planetary satellites and the 
terrestrial planets, in which the dissipation decreases when the frequency increases (Efroimsky-Lainey 
regime), and the other, including giant planets, hot Jupiters and stars, following Darwin's regime in 
which dissipation increases with the frequency. 

At last, the Appendix discusses some accuracy problems involved in tidal evolution theories. The 
warning included in FRH stating that long series expansions by themselves are not improvements of a 
physical theory has been several times misunderstood and criticized (e.g. Leconte et al. 2011). However, 
we repeat it. Very-high-order expansions did create, in the past, the feeling that the results of standard 
Darwin theory were exact notwithstanding the fact that we cannot yet say that we understand the 
Physics ruling bodily tides. Indeed, we guess that the results of the standard Darwin theory for the 
pseudo-synchronous stationary rotation of close-in companions, independent of the viscosity of the 
body, cannot be correct, even when written using endless series. 

On the negative side, the main question is that the rheophysical theory failed to give the actual 
shape of the bodily tide observed in the Earth. This failure forced us to admit the existence of a 
superposed elastic tide, which affects the shape of the observed tide but which is torque free and thus 
does not change the results of the rheophysical theory in what concerns the rotation of the bodies and 
the average energy dissipation. This superposition of the creeping tide and an elastic tide is a question 
that certainly needs to be stated in terms of Physics. 

In addition, we can devise many other points in which the theory needs to be improved. We list a 
few of them: We have considered only homogeneous bodies; we have considered a non-rotating model 
for the equilibrium tide (a Jeans spheroid); we have considered only the planar (two-dimensional) 
problem; we have disregarded inertia and the non-linearity of actual creep laws; we have assumed tidal 
deformations small enough so as to allow a superposition principle to be used in the calculation of 
the potential spanned by the tidally deformed body, etc. In that sense, the presented theory is yet 
a proposed model. Each of the cited points deserves now to be taken into consideration, and shall 
be taken into consideration if we want to apply the theory to bodies as intricate as some planetary 
satellites. One positive point to be mentioned is that the theory opens the way for the construction of 
very complex models and to adopt laws more complete than the Newtonian creep law. Once the basic 
equations are given, they can be solved numerically. This means that we may adopt physical models 
as complex as necessary, no matter if their equations can be solved analytically or not. 

At last, we have to emphasize that the approach developed in this paper is a new and complete 
theory of the bodily tide problem, whose results derive from only one physical law: the Newtonian creep. 
There are no ad hoc lags plugged by hand as in the standard theories. The constants (Tj appearing in 
the arguments of the trigonometric functions, which can be physically interpreted as delays, are well 
determined parameters of the exact solution of simple ordinary differential equations. The results show 
some discrepancies with respect to the classical theories, but also some coincidences. The way in which 
creep and elastic tides combine to give rise to geodetic lags increasing when the frequency decreases, 
for hard bodies is appealing and may serve as a justification to the modern theories of Efroimsky and 
collaborators. In what concerns the discrepancies, it is important to know if they refer to observable 
phenomena or not. When observational data do not exist, it is impossible to know whether a proposed 
theory is or not correct. 
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A Higher-order approximations 

The eccentricity functions introduced in section [3] are the Cayley expansions for the solutions of the Keplerian 
motion (see Cayley, 1861; da Silva Fernandes, 1996) corresponding to the entries (r/a) _3 [sin | cos]2/ in Cayley's 
tables. 

The approximations adopted in Sec. [3] may be enough for most purposes. They are certainly enough for 
low eccentricities. For high eccentricities, however, a different approach is necessary. Cayley expansions include 
in their derivation the Taylor expansion of the solution of Kepler's equation, whose convergence radius is 
e* = 0.6627434... (see Wintner, 1941). Some simple comparisons allow us to see that for eccentricities of the 
order of 0.5, or even less, the results are not accurate. 



Table 2 Comparison of the calculated values of i?2,fe(e) for e = 0.4 
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Since the equations appearing in this pa per are written in terms of one of the families of Cayley expansions, 
the power series may be substituted by eq. \W\ in the form 

E 2 k (e) = ; / - cos (2v + (k - 2)1) dv. (63) 

2-rry/l - e 2 Jo r 

We may note that the function under the integral involves both the mean (£) and the true (v) anomalies, but 
can be easily handled numerically. The only difference with eq. ([16[) is the preference in having the integration 
done over the true anomaly instead of being done over the mean anomaly (via the classical expression r 2 dv — 
a 2 y/l — e 2 dt) to avoid having to solve the Kepler equation, which is operationally expensive. 

The integral solves the problem of the accuracy in the determination of the Cayley coefficients, but not that 
of the Fourier series giving the variation of the elements. To obtain approximations valid in high eccentricities, 
it is necessary to use expansions not affected by the singularities of the Keplerian motion (see Ferraz-Mello and 
Sato, 1989). However, one preliminary guess can be done comparing different approximations. Fig. [8] shows the 
value of \e\ for close-in companions in pseudo-synch rono us stationary rotation, in arbitrary units, considering 
the sum of the eccentricity-dependent terms of Eqn. [49] truncated at TV = 7, 20, 50 and 100, respectively. The 
abrupt changes of the curves with respect to the next one show a limiting eccentricity in each case. For instance 
the series truncated at TV = 7 (the usual cutoff of Cayley series) give wrong results for eccentricities above 
0.4. It is important to stress that the calculations showed in fig. [5J concerning the tidal evolution predicted by 
the standard Darwin theory are not affected by truncation effects since for that sake we have used Mignard's 
formulation whose basic equations are given in closed form (Mignard, 1979) without resorting to any expansions. 
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